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Foreward 


This  report  is  Chapter  XII  of' the  twelve  in  a  forthcoming  research 
monograph  on  the  mathematical  theory  of  laminar  combustion.  Chapters  I-IV 
originally  appeared  as  Technical  Reports  Nos.  77,  80,  82  &  85;  these  were 
later  extensively  revised  and  then  issued  as  Technical  Summary  Reports 
No's  l803,  l8l8,  l8l9  &  1888  of  the  Mathematics  Research  Center,  University 
of  Wisconsin-Madison.  References  to  I-IV  mean  the  MRC  reports. 
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Chapter  XII 
Ignition  St  Explosion 

1.  Synopsis. 

The  phenomenon  of  ignition  has  already  "been  encountered  in  at  least 
three  contexts:  "burning  of  a  linear  condensate  (Ch.  IV) ,  spherical  diffusion 
flame  (Ch.  VI)  and.  spherical  premixed  flame  (Ch.  VII).  From  the  S-shaped 
response  curves  determined  by  steady-state  analyses  it  was  argued  that  the 
burning  rate  would  jump  from  a  weak,  almost  extinguished,  level  to  a  vigorous 
one  as  the  pressure  (and  hence  the  Damktfhler  number)  is  increased  through 
some  critical  value  (corresponding  to  the  lower  bend  of  S).  An  intrinsically 
unsteady  phenomenon,  appropriately  called  ignition  in  these  contexts,  was  . 
thereby  inferred  from  existence  results  for  the  steady  state. 

Such  analyses  have  a  fairly  long  history  in  combustion,  the  simplest 
example  being  the  .thermal  theory  of  spontaneous  combustion  identified 
with  the  name  of  its  originator,  Frank-Kamenetskii  (l 969).  Sec.  2  will  give 
a  mathematical  version  of  that  theory,  which  shows  an  early  appreciation  of 
activation-energy  asymptotics  (though  not  in  the  formal  sense  of  the  present 
monograph). 

Useful  though  steady-state  analyses  undoubtedly  are,  such  phenomena 
demand  unsteady  descriptions;  and  that  is  the  main  goal  of  the  present 
chapter.  The  three  contexts  already  encountered  are,  however,  too 
complicated  for  our  purposes  and  we  turn  instead  to  a  simpler  problem 
containing  the  essential  feature,  namely  the  evolution  of  a  deflagration 
wave  from  an  unburnt  state. 
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Consider  a  thermally  insulted  enclosure  containing  a  mixture  at 
sufficiently  low  temperature  for  the  reaction  though  present,  to  be  very 
weak.  Since  the  heat  released  cannot  escape  the  temperature  must  rise, 
albeit  very  slowly.  After  a  long  time  the  heat  generated  by  the  initially 
weak  reaction  will  have  raised  the  temperature  of  the  mixture  enough  to 
excite  rapid  reaction  somewhere  and,  if  the  combustion  field  is  inhomogeneous, 

« 

a  deflagration  wave  win  sweep  across  it.  To  distinguish  this  type  of 
rapid  reaction  from  that  caused  by  external  agencies,  the  term  thermal 
explosion  is  used.  (Thermal  refers  to  the  temperature  sensitivity  of  the 
reaction  rate.  Another  type  of  explosion  is  due  to  a  rapid  chain  reaction 
being  initiated  when  conditions  are  such  that  a  .few  necessary  radicals  are 
produced  from  one  of  the  reactants.) 

A  spatially  homogeneous  explosion  is  described  by  ordinary  differential 
equations.  Its  study  in  Sec.  3  provides  insight  into  the  mathematical  nature 
of  the  rapid  transients  characteristic  of  the  explosion  process.  That 
prepares  the  way  for  a  discussion  of  the  spatially  inhomogeneous  problem 
in  Sec.  4. 

A  distinctive  feature  of  thermal  explosion  is  the  very  rapid  temperature 
increase,  known  as  thermal  runaway,  that  occurs  at  the  end  of  a  well-defined 
induction  period.  The  mathematical  manifestation  is  unboundedness  in  the 
solution  of  an  equation  for  small  perturbations  of  the  initial  state. 
(Evolution  equations  that  generate  unboundedness  solutions  after  a  finite 
time  have  recently  attracted  the  attention  of  mathematical  analysts.)  The 
subject  of  tbsa  thermal  ignition  provides  other  examples,  one  of  the  simplest 
being  a  half  space  filled  with  a  combustible  mixture,  at  the  boundary  of 
which  a  positive  heat  flux  is  applied.  After  a  certain  time,  ignition 


occurs  at  the  "boundary  and  a  plane  deflagration  wave  propagates  into  the 
interior,  consuming. the  mixture  in  its  passage.  The  early  stages  of  the 
process  are  described  in  Sec.  5. 

The  problems  in  Secs.  2-5  are  all  concerned  with  reactants  that  are 
premixed;  a  whole  range  of  such  problems,  to  which  activation-energy 
asymptotics  could  be  profitably  applied,  is  discussed  by  Merzhanov  & 

Averson  (1971).  By  contrast  Sec.  6  treats  the  so-called  Marble-Adamson 
(195*0  problem,  in  which  a  fuel  and  oxidant  in  separate  half-spaces  are 
brought  into  contact  at  the  initial  instant  so  as  to  interdiffuse  and  burn. 
There  are  many  points  of  similarity  with  the  spatially  inhomogeneous  thermal 
explosion,  and  a  reasonably  complete  description  of  the  combustion  process 
is  possible. 


Spontaneous  Combustion. 


Exothermic  chemical-reaction  is  a  well-known  hazard  of  certain  bulk 
materials  in  storage  or  transit.  Joseph  Conrad  in  his  story  "Youth" 
graphically  describes  the  peril.  Take  one  ship  with  a  cargo  of  coal 
dampened  down  because  of  an  earlier  leak  (now  fixed);  embark  on  a. slow 
voyage  to  Bangkok;  and  before  the  China  Sea  is  reached  the  legend  "Do  or 
Die"  proves  ominous  rather  than  brave. 

Any  moist  volume  of  organic  material  (e.g.  wood  chips,  wool,  corn 
cobs,  compost)  will,  because  of  chemical  reaction  initiated  by  micro¬ 
organisms,  get  hot.  (The  reduction  of  spontaneous  reactions  in  fruits 
and  vegetables  is  a  major  concern  of  post-harvest  physiologists.)  Certain 
inorganic  material  such  as  ammonium  nitrate,  lead  azide  or  a  hydrogen  sulphide/ 
oxygen  mixture  will  also  heat  up  because  of  chemical  decomposition. 


If  the  reaction  were  independent  of  temperature  there  would  he  a 

maximum,  depending  on  the  dimensions,  to  which  the  temperature  of  a  volume 

of  material  would  rise  after  a  long  time.  (That  maximum  holds  steady  over 

an  even  longer  period,  but  there  is  an  eventual  decay  due  to  reactant 

depletion.)  In  the  quasi-steady  state  the  rate  of  generation  of  heat, 

which  is  proportional  to  a  for  a  volume  of  characteristic  dimension  a, 

would  be  balanced  by  loss  through  the  surface ,  which  is  proportional  to 
2 

a  with  a  factor  of  proportionality  an  increasing  function  of  temperature 
there.  Clearly  the  larger  the  volume'  in  equilibrium  the  hotter  it  would  be. 
However,  the  reaction  does  depend  on  temperature,  proceeding  much  more 
vigorously  at  higher  temperatures ,  so  that  there  will  be  a  critical  value  of 
a  beyond  which  the  efflux  at  the  surface  is  insufficient  to  dispose  of  the 
heat  generated  and  unsteady  combustion  will  occur.  Our  first  task  then 
is  to  look  for  steady  states. 

Any  steady  state  is  described  in  equation  (I.56)  with  8/3t  =  y  =  0,  i.e. 
(1)  '  V2T  *  -ft  where  ft  =  De"6^T, 


the  latter  being  appropriate  when  the  length  a  is  used  in  non-dimensior.aliza 

tion.  For  solids,  the  density  has  a  constant  value,  which  may  be  designated 

Pc  so  as  to  replace  p  by  1  in  the  definition  (1.5 9)  of  D;  pressure  does 

not  arise,  which  means  that  p  should  be  replaced  by  p  QR/c  m  in  that 

c  c  p 

definition  (thus  eliminating  the  gas  constant  R  from  it) ;  and  a  is 
usually  taken  to  be  zero.  The  only  new  assumption  is  that  reactant  depletion 
may  be  neglected  (the  constant  mass  fractions  in  ft  being  absorbed  into 
T>)  because  the  steady  state  is  close  to  the  initial  state,  a  condition  ensured 
a  posteriori  when  0  is  large.  For  gases,  of  primary  interest  to  us, 
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there  is  the  additional  assumption  of  small  heat  release  so  that  the  constant 
density  approximation  can  be  made  to  ensure  that  the  mixture  stays  effectively 
at  rest. 

For  an  initially  uniform  temperature  Tq  such  as  we  have  so  far 
envisaged,  this  purely  thermal  problem  can  be  treated  in  the  limit  0  » 

by  setting 


(2)  T  =  TQ  +  e"1^  +A 

To  »m£  order  we  find 


(3) 


V2$  a  r5e*, 


where 

(4)  *  *  Tx/T2  and  6  *  D0e  °/T2. 

Note  that  the  steady  state  being  sought  is  close  to  the  initial  state, 
so  that  neglect  of  reactant  depletion  is  justified.  However,  for  natural 
organic  materials  such  as  wood  chips  the  activation  energy  is  not  large 
and  depletion  is  not  small;  in  either  circumstance  the  present  analysis 
is  not  appropriate. 

The  approximation  (2)  is  familiar  from  our  earlier  ignition  studies 
(Secs.  VI. 4  &  VII. 6).  In  the  present  context  it  is  due  to  Frank-Kamenetskii , 
constituting  an  early  example  of  activation-energy  asymptotics.  References 
and  a  discussion  of  the  solutions  of  equation  (3)  are  given  in  his  book 
(1969).  On  the  other  hand,  linearization  of  the  exponent  0/T  in  the 
Arrhenius  factor  has  sometimes  been  adopted  as  a  model  not  rationally 
consistent  (in  an  asymptotic  sense)  with  the  underlying  equation  (l)  for  t 
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the  specific  problem  being  examined.  The  plane  form  of  equation  (3)  has 
a  long  history  (Stuart  (1967),  going  back  at  least  to  Liouville  (1853),  to 
whom  it  is  sometimes  attributed.  The  general  solution  is  known,  but  there 
appears  to  be  no  way  of  fitting  it  to  boundary  conditions  of  interest.  As 
a  consequence,  solutions  of  complete  boundary-value  problems  have  only 
been  constructed  in  a  few  special  cases. 

Amongst  these  are  problems  depending  on  x  alone.  Indeed  the  correspond¬ 
ing  problems  for  the  exact  equation  (l)  can  then  be  solved,  providing  useful 
checks  for  the  present  asymptotics  [cf.  Shouman  fc  Donaldson  (1975)]*  A 
particularly  simple  example  is  an  infinite  slab  |  oc  j  <_ 1  of  combustible 

material  whose  surfaces  are  held  at  the  temperature  T  ,  so  that 

o 

(5)  if>  =  0  for  x  =  ±1. 

Such  a  surface  condition  corresponds  to  the  limit  of  an  efflux  k(?  -  ?Q) 
per  unit  area  as  k  ->■  »,  and  it  is  the  *  ility  of  the  material  in  the  slab 
to  conduct  that  heat  to  the  surfaces  whxth  determines  whether  a  steady  state 
exists  or  not.  By  symmetry  <£  will  attain  its  maximum,  say,  at 

x  =  0-  and  in  terms  of  it  we  have 

(6)  $  =  2  in  ( e  m  sech^e^/2^'“x)] . 

(Cf.  Fig.  2a  for  a  graph  of  this  function.)  The  boundary  conditions 
(5)  therefore  require 


a  relation  fixing  the  maximum  temperature .in  terms  of  6.  For  a  given 
material  6  can  be  varied  through  a,  the  semi -thickness  of  the  slab, 


.  tii.- ^  ^  ,  -XSUir-. 
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or  Tq,  the  temperature  at  its  surfaces  [see  equation  (kb)  and  the  definition 
(1.60)  of  D]. 

Fig.  1  shows  the  response  (7),  which  has  the  typical  shape  of  an 
ignition  curve.  There  is  a  critical  value 

(8)  5  =  .878. 

c 

above  which  no  solution  exists  and  below  which  there  are  two  solutions .  A 
value  of  5  greater  than  6c  is  said  to  be  supercritical  while  a  smaller 
value  is  said  to  be  subcritical.  The  implication  is  that  thick  hot  slabs, 
being  supercritical,  will  burst  into  flames  but  thin  cold  slabs,  being 
supercritical,  will  not;  and  this  general  picture  has  been  confirmed  experi¬ 
mentally  for  a  number  of  substances. 

The  symmetry  of  the  boundary  conditions  makes  the  solution  atypical  however, 
as  we  shall  see  from  another  example  in  which  the  surfaces  x  =  +1  are  held 
at  temperatures  TQ  and  aTQ,  respectively,  where  0  <  a  <  1.  Now  the 
reaction  is  confined  to  a  neighborhood  of  the  hot  boundary  in  the  limit 
6  ■+■  »  instead  of  being  spread  across  the  slab,  so  that  conduction  of 
its  heat  to-  that  surface  is  much  easier  and  is  much  larger. 

The  initial  temperature  in  the  slab,  due  entirely  to  conduction,  is 

(9)  T  s  [l  +  a  +  (1  -  a)x]TQ/2 

and  we  seek  a  steady  state  within  O(0_1).  For  such  a  state  the  0(l) 
temperature  has  its  maximum  at  the  hot  boundary,  so  that  reaction  only 
occurs  in  a  2one  near  there.  The  expansion  (2)  is  still  valid,  but 
its  variable  is  now  • 


(10) 


5  -  8(1  -  x) 


The  structure  of  the  reaction  zone  is  found  to  "be  governed  by 


(11) 


d2<!>/d52  =  -  de*  with  5  =  0_2d 


and 


(12) 


d<^/d£  =  (l  -  o)/2Tq  +  b(l)  as  -»,$  =  0  for  5=0, 


where  the  first  boundary  condition  comes  from  matching  with  the  temperature 
(9)  outside  the  zone.  The  general  solution  of  equation  (ll)  is 


(13) 


*0/2 

$  =  2  in  sech 


(5-  50)l} 


where  5Q  is  the  location  of  the  maximum  and  to  satisfy  the  conditions 


(12)  23  must  sei 


(lU) 


5  =  (l-g)  e  *°  and  5  =  ±(— |r-)1/2  sech"1{e 


Since  <J>  vanishes  at  the  surface,  ranges  from  0  to  »  (otherwise 


50  is  complex) ,  so  that  6  cannot  exceed 

(15) 


6C  =  (1  -  o)2/8T ®. 


The  +  signs  on  5Q  provide  two  solutions  when  6  <dc  (see  Fig.  2a);  for 
the  plus  sign,  5Q  lies  outside  the  slab  and  the  maximum  t.  in  the  slab 
is  zero  (its  surface  value)  and  not  If  denotes  this  maximum,  then 


"m 


(16) 


rm 


=  f  ° 

Vo 


♦a  for  5q  <  0, 


for  50  >  0 


and  the  relation  (lH)  leads  to  the  response  shown  in  Fig.  2b. 


„,.  ,  i r ■**■■  ■■«  >.'_■  "  ■«■—«—»'■■  ~ .  J  / -  11  ^ 
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Note  that  6=0  6  implies  that  the  D  corresponding  to  6^  is  much 

larger  than  that  corresponding  to  6q.  We  conclude  that  raising  the 

temperature  of  the  cold  boundary  to  that  of  the  hot  will  greatly  increase 

the  risk  of  explosion.  In  fact  the  formula  (15)  suggests  that  the  risk 

increases  sharply  as  soon  as  o  is  within  0(  6**1)  of  1. 

Two  further  points  should  be  noted.  The  second  problem  illustrates 

ignition  rather  than  explosion.  With  an  ambient  temperature  ciTq,  an 

external  agency  is  required  to  raise  the  surface  temperature  at  x  =  1 

to  T^;  and  if  that  is  ultimately  done  slowly  enough  the  conduction  profile 

(9)  will  be  established.  Also  we  note  that  nothing  has  been  missed  in 

either  problem  by  looking  for  steady  states  within  0(e”'1')  of  the  initial 

state.  Any  steady  state  further  away  would,  correspond  to  a  D  of  exponential 

order  e/T  ,  where  the  maximum  temperature  T  is  greater  than,  Tft.  It 
m  m  o 

would  therefore  correspond  to  the  origin  in  Figs.  1  and  2b. 

For  a  review  of  the  subject,  including  other  geometries  and  boundary 
conditions,  see  Gray  &  Lee  (1967).  They  do  not  use  formal  asymptotics, 
and  indeed  there  is  room  for  a  comprehensive  treatment  doing  so. 

3.  Adiabatic  Explosion. 

To  describe  an  unsteady  phase  of  the  combustion,  equation  (l)  is 
replaced  by 

(17)  9T/3t  -  V2T  =  De-0^ 

so  that  the  small-disturbance  equation  (3)  becomes 

(18)  3 $/  at  -  V2$  =  6e^. 
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Th-  problem  is  completed  by  the  initial  condition 

(19)  <j>  =  0  for  t  =  0 
and  suitable  boundary  conditions. 

The  non-existence  of  a  steady-state  solution  for  the  slab  with  equally 
hot  surfaces  and  supercritical  6  implies  that  $  must  increase  without 
bound,  i.e.  that  T  goes  further  than  0(e“1)  from  TQ.  Indeed  this  happens 
at  a  finite  time,  as  we  shall  now  show  for  the  spatically  homogeneous  problem 
of  a  material  thermally  insulated  at  its  boundary  surface,  when  equation  (.17) 
reduces  to 

(20)  dT/dt  =  De“9^T. 

Our  object  in  the  present  section  is  to  gain  insight  into  thermal  explosions 
generally  by  examining  those  that  are  spatially  homogeneous. 

The  corresponding  small-disturbance  equation  is 

(21)  d<f>/dt  =  6e^, 
which  has  the  solution 

(22)  <j>  =  -£n(l  -x) 
satisfying  the  initial  condition  (l9)j  here 

(23)  x  =  fit. 

At  the  finite  induction  time  6  \  the  perturbation  <}>  becomes  unbounded 
and  we  speak  of  thermal  runaway..  Note  that  the  phenomenon  occurs  for  all 
values  of  6  in  the  spatially  homogeneous  problem  because  the  heat  generated 
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by  the  reacts^  cannot  escape,  i.e.  conditions  are  adiabatic.  On  the 
other  hand,  if  an  artificial  heat-loss  term  is  introduced  then  subcritical 
behavior  can  be  restored,  see  Kassoy  &  Poland  (1975). 

Clarification  of  the  runaway  process  requires  a  return  to  equation  (20). 
However,  beyond  thermal  runaway  a  purely  thermal  model  is  of  little  practical 
significance  since  it  predicts  an  unbounded  increase  in  T,  whereas  the 

rise  ceases  when  the  reactant'  is  consumed.  Under  the  assumption  that  a 

* 

single  (decomposing)  reactant  is  involved,  its  depletion  can  be  incorporated 
by  reinstating  the  mass  fraction  in  equation  (20),  i.e.  writing 

(24)  dT/dt  =  D(Y/Yo)e"0//,r, 

and  adding 

(25)  dY/dt  =-D(Y/Yo)e"0/T, 

where  YQ  is  the  initial  fraction  of  reactant  in  the  material.  Clearly 
the  Shvab-Zeldovich  relation 

(26)  Y  +  T  =  YQ  +  .T0  =  Kq  (say) 

holds.  No  change  in  our  analysis  of  the  induction  phase  (Sec.  2)  is  necessary 
because  Y  then  stays  within  0(0  1)  of  Y  and  hence  may  be  replaced  by 
it  in  equation  (24).  We  are  therefore  led  to  the  problem 

(27)  dT/dt  =  (TqS/YqQ) (Hq  -  T)exp(0/To  -  o/T) ,  T  =  TQ  for  t  *  0. 

Clearly  T  is  monotonically  increasing  but  the  factor  (hq  -  T-),  corresponding 
to  reactant  depletion',  ensures  that  it  does  not  go  beyond  in  fact  it 

rapidly  approaches  Hq  (exponentially  in  time). 
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’  Kassoy  (1975)  hag  discussed  the  asymptotics  of  this  problem  directly; 

here  we  shall  verify  the  salient  features  of  his  results  from  the  exact 
solution.  With  the  confidence  so  gained  we  shall  attack  spatially 
inhomogeneous  problems  in  Sec.  4,  for  which  there  are  ho  exact  solutions 
against  which  to  check  the  asymptotics.  The  exact  solution  here  is  expressed 
in  terms  of  the  exponential-integral  function 

(28)  Ki(x)  =  |  % 

.  —CO 

where  the  principal- value  sigh  is  unnecessary  for  x  ^  0.  With  the  time 
variable  (23)  we  find 

0/t  q/h 

(29)  Tq  e  °(xe-T  /YqQ  =  e  0  Ei(e/T  -  0/Hq)  -  Ei(e/T), 

>  'where- 

_Q/T  0/h 

(30)  xe  =  (YQ0e  °/Tq)  [e  Sito/T^e/H^-rEite/^)]. 

Since 

(31)  Ei(x)  =  £n|x|  +  0(l)  as  x  ->  0 

the  formula  confirms  that  T  approaches  Hq  exponentially  as 
0~1exp(0/To-0/HQ)  t  -*». 

Three  phases  can  be  distinguished  when  0  is  large:  induction, 
explosion  and  relaxation.  If  T  is  bounded  away  from  H^,  the  arguments 
of  both  Ei-functions  in  the  solution  (29)  become  large  so  that,  since 

,  (32)  Ei(x)  =  x  ■*’eX[l  +  x  ■*'  +  0(x  ^)]  as  x  ->  +«, 

we  find 
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(33)  re-T  =  ryoT2/T|(Ho-T)]exp(0/T-0/To)[l+O(0"1)]> 

which  reveals  two  of  the  phases.  In  the  first  phase  T  is  still  within 
0( s”"1)  of  and  the  right-hand  side  is  0(l),  corresponding  to  induction. 

In  fact,  the  result  (22)  is  recovered  when 

(34)  xe  =  1  +  TQ(2  +  To/yo)0_1  +  O(0“2) 

is  noted.  In  the  second  phase  T  is  hounded  away  from  (as  well  as 

from  Hq)  and  the  right-hand  side  is  exponentially  small,  showing  that 
most  of  the  change  in T occurs  within  an  exponentially  small  interval  just 
before  the  time  Tg*  In  that  sense  Tg  is  the  time  at  which  the  explosion 
takes  place. 

The  third,  relaxation  phase  corresponds  to  values  of  T  within 
/  -1\ 

0(0  )  of  Hq;  then  the  argument  of  the  first  Ei-fiinction  in  the  solution 

(*?)  is  0(l)  and  we  have 

(35)  t6-t  *  (Yo0/T2).exp(0/Ho-0/To)Ei(0/T  -  0/HQ)[l  +  0(e'1)].' 

The  final  consumption  of  the  reactnat  also  takes  place  within  an  exponentially 
small  interval  which,  since  the  Ei-function  changes  sign  at  some  positive 
values  of  its  argument,  spans  x  . 

We  now  turn  to  the  question  of  deriving  these  asymptotic  results  directly 
from  the  governing  equation  (24).  Consider  first  the  induction  phase  and 
set 

(3 6)  t  =  tq  +  e_1T1(t)  +  0”2T2(t)  +  o(e"3). 

Then  we  find 


as  in  the  earlier  result  (22),  and 

(38)  T2  =  Tq  Un2(l  -  t)  -  (2  +  T0/YQ)[to(l  -  x)  +  t/(l  -  x)]}. 

The  logarithm  in  T^  suggests  that  breakdown  occurs  when  t  is 
exponentially  close  to  1,  while  the  pole  (l  -  x)~^  in  Tg .  changes  that 
to  algebraically  close;  and  correctly  so  since  the  explosion  time  (34)  is 
6(0_1)  away  from  1.  To  restore  exponential  closeness  it  is  only  necessary 
to  absorb  the  pole  into  Tn  by  writing 

(39)  Tx  =  -£n[l  +  Tq(2  +  T0/Y0)o"1  -  x], 

which  amounts  to  introducing  the  0(8*"1)  approximation  to  xg.  The  same 
end  can  be  achieved  formally  by  considering  both  T  and  x  to  be  functions 
of  a  new  time  variable  x*.  If 

(40)  x  =  x“+  0  1T1(x')  +■  6(0  2) 
with 

(41)  x-^Oj  =  0  and  x3(l)  =  TQ(2  +  Y  0/Tq) 

(  x^  otherwise  being  arbitrary),  then  the  explosion  is  fixed  at  x  =  1 
and  no  pole  occurs  in  T2(x).  In  short,  the  time  to  explosion  has  been 
determined  asymptotically  by  considering  only  the  induction  phase;  it  is 
reasonable  to  suppose  that  the  same  procedure  will  work  in  other  problems. 

The  description  (33)  of  the  explosion  phase  can  also  be  derived  directly 
from  equation  (24)  but  we  shall  not  go  into  details.  The  most  important 


feature  is  the  time  scale 


(42)  a  =  -e_1£n(t-  -  t), 

e 

which  plans  a  key  role  in  the  study  of  the  inhomogeneous  problem,  as 
we  shall  see  in  Sec.  4.  As  a  ranges  from  0  to  Y^/T^H^  the  temperature 

(43)  T  =  TQ/(l  -  oTq)  +  0(1) 

ranges  from  TQ  to  H^,  to  leading  order. 

For  the  final,  relaxation  phase  we  write 

(44)  t  «  hq  +  e"1^**)  +  a(e"2) 

in  equation  (24),  which  shows  that  the  appropriate  time  variable  is 

(45)  V  =  (t  -  t  )/-$ 

e 

where 

(46)  C  *  (Yo6/T2)exp(-0Yo/HoTo). 

Then 

(47)  d<j>/d»  =  with  $  =  T-jH2, 
the  general  solution  of  which  is 

(48)  u  -  ^  =  Ei(l)  -  Ei(-<j>); 

here  ^  is  just  the  value  of  a  when  $  =  -1.  The  requirement  that  <j>  -M3 
as  «  -»•  +®  (i.e.  all  reactant  is  consumed)  is  satisfied  automatically,  as 
is  that  of  matching  with  the  expansion  (43).  For  the  latter,  note  the 


result 
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(1*9)  »  =  e  ^/4» 

of  introducing  <*  into  the  expansion  (33)  and  then  keeping  it  fixed  under 
expansion.  Comparison  with  the  solution  (H8)  shows  that  the  integration 
constant  is  left  undetermined  and  in  fact  cannot  he  determined  by 

asymptotic  development:  however  many  terms  are  retained  in  the  expansion 
(33)  it  will  never  recapture  the  Ei-functions  of  '..he  exact  solution,  whereas 
those  functions  determine  <j>.  It  is  therefore  fortunate  that  results 
analogous  to 

(50)  W1=-Ei(l), 

here  obtained  from  the  exact  solution,  will  not  be  needed  in  the  sequel. 

4.  Explosion  with  Heat  Loss 

Kassoy  ( 197T )  has  incorporated  non-adiabaticity  into  the  spatially 
homogeneous  problem  by  inserting  a  distributed  heat-loss  term  proportional 
to  T  -  Tq  into  equation  ( 2k ).  Apart  from  the  lack  of  a  Shvab-Zeldovich 
relation,  the  new  features  are  the  inability  of  the  reaction  to  produce 
the  adiabatic  temperature  HQ  and  the  eventual  decay  to  TQ.  To  obtain 
the  essentially  new  phenomena  of  hot-spot  formation  and  subsequent 
flash-through  it  is  necessary  to  consider  a  more  realistic  problem 
involving  spatial  inhomogeneity,  and  that  can  be  done  without  sacrificing 
the  Shvab-Zeldovich  relation. 

Consider  once  more  the  infinite  slab  with  equally  hot  surfaces. 

The  equations  ( 2h  ,2!j)  are  now  replaced  by 


(51) 


3T/3t  -  3  T/3x2  =  -Oy/3t  -^Vy/ax2)  =  p(Y/Y  )e 


-e/T 


i 


cf.  equation  (1-7).  For  a  solid,  where  the  reactant  is  unable  to  diffuse, 
the  Lewis  number  is  infinite.  Nevertheless  we  shall  take£  =1  so  as 


to  have  the  simplification  of  a  Shvab-Zeldovich  integral  (26)7?  Rather  than 
maintaining  the  reactant  there  it  would  be  more  practical  to  prevent  its 
flux.  We  shall  see,  however,  that  the  boundary  condition  on  Y  does  not 
influence  the  solution  until  after  the  main  events  have  taken  place.  The 
problem  we  shall  treat  is  therefore 


3T/3t- -a 2t/3x2  =  (T26/Yo0)(Ho-T)exp(0/To  -  q/ T) 


T  =  T.  for  t  =  0  and  for  x  =  ±1. 


x 


By  symmetry  the  condition  at  x  =  -1  can  be  replaced  by 


(5M 


3T/3x  =0  for  x  =  0, 


with  only  the  half-slab  0  x  <_  1'  considered. 

In  addition  to  the  induction,  explosion  and  relaxation  phases  found 
in  the  spatially  homogeneous  problem  (27)  there  are  now  transition  and 
propagation  phases.  These  arise  because  the  explosion  is  here  confined  to 
a  single  location  and  then  propagates  elsewhere,  involving  a  transition 
from  the  formation  of  a  hot  spot  to  the  flash-through  of  a  deflagration 
wave.  These  two  new  phases  postpone  the  relaxation  phase,  which  is  the 
only  feature  affected  by  the  boundary  condition  on  Y. 

As  for  the  homogeneous  problem,  the  induction  phase  is  described  by 
an  expansion  of  the  form  (36),  but  nov  the  coefficient  functions  depend 


on  x  also  and  T,  =  T satisfies 
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(55)  3<f>/3t  -  32<J>/3x2  =  de^, 

(56)  3<j>/3x  =  0  for  x  =  0,  <f>  =  0  for  x  =  1  and  for  t  =  0. 

If  6  is  greater  than  H ,  as  we  shall  suppose,  the  solution  cannot  tend 
to  a  steady  limit.  Since  the  solution  of  the  corresponding  homogeneous 
problem  becomes  unbounded  in  a  finite  time,  it  is  not  surprising  to 
find  that  the  same  is  true  here  (cf.  Payne  1975  and  Ayeni  1978).  Thermal 
runaway  can  no  longer  occur  everywhere  when  the  surfaces  are  held  at  T^  but 
is  localized  at  the  symmetry  plane  x  «  0.  The  runaway  time  d  \  ( 6)  is  no 
longer  fixed  on  the  x-scale,  but  is  a  function  of  6  that  was  first 
determined  (numerically)  by  Poland  (1979),  see  Kassoy  &  Poland  (1980a) , 
and  then  independently  by  Kapila  (198o).  As  6+  6C  the  time  to  runaway  becomes 
indefinitely  long. 

While  <f>  cannot  be  calculated  analytically,  its  form  as  x  ->•  xg  can  be 
determined,  a  step  that  is  crucial  for  the  subsequent  development.  Backward 
integration  of  a  parabolic  equation  away  from  a  singularity  (in  this  case 
at  t  =  6  \  )  is  no  novelty;  more  than  thirty  years  ago  Goldstein  (19^8) 
encountered  such  a  question  in  his  discussion  of  Prandtl's  boundary-layer 
equations  near  a  stagnation  point.  He  showed  that  the  similarity  variable, 
here 

(57)  n  =  6l/2x/(xe  -  x)l/2, 

plays  a  central  role  near  the  singularity.  In  terms  of  q,x  our  equation 
becomes 

(58)  3<(>/3x  -  (xe  -  x)  32<f>/3n^  “  (n/2)3<f>/3ri3  =  e*. 
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In  view  of  the  result  (.3*0  for  the  homogeneous  problem,  a  uniform 
asymptotic  appxoximation 

(59)  4>  =  -  £n(x  -  t)  +  <Js  (n)  +  «(l)  as  x  +  t 

c  U  v 

is  sought  for  n  bounded.  We  find 

*0 

.(6o)  -  (n/2)^  +  e  =1  where  <^(0)  =  0 

according  to  the  symmetry  condition  (5*0;  the  other  boundary  conditions 
cannot  be  applied.  But  to  determine  <j>g  another  boundary  condition  is 
needed  and  that  comes  from  matching  with  the  induction  solution.  Since  ^ 
cannot  be  exponentially  large,  it  has  the  asymptotic  form 

(61)  <j>Q  =  -2fcnn  +  A  +  oil)  as  n  -►  *  > 

which  provides  the  needed  condition  once  the  constant  A  has  been 
determined  by  matching.  To  effect  the  matching  we  write  the  explicit 
terms  (6l)  as  functions  of  x  and  x,  and  expand  in  xg  -  x  to  obtain 

(62)  <(>=  -2  in  x  -  in  6  +  A  +  o(l) 

as  the  behavior  of  the  induction  solution  as  x  f  xg  with  x  bounded  away 
from  zero.  The  numerics  mentioned  above  do  indeed  exhibit  such  a  behavior, 
so  that  A(d)  can  be  determined  (Poland  1979,  see  Kassoy  &  Poland  1980a; 
Kapila  1980).  The  procedure  fails  if  6  is  subcritical  because  the 
numerical  integration  does  not  yield  the  required  form  (62)  but  tends  to  a 
steady  state  instead. 

The  expansion  (59)  exhibits  a  focussing  effect  in  that  constant  values 
of  moves  towards  x  =  .0  as  x  t  igl  i.e.  the  peak  of  the  temperature 


distribution  for  t  constant  becomes  sharper  as  thermal  runaway 
approached.  This  feature. persists  into  the  explosion  phase,  as  ve  shall 
see  shortly,  so  that  a  well-defined  hot  spot  is  formed.  A  conibinat.se 
asymptotics  and  numerics  therefore  provides  a  complete  picture  for  t  < 
Induction  is  immediately  followed  by  explosion,  discussed  by  Kassoy  & 
Poland  (1980b)  and  by  Kapila  (1930),  for  which  the  time  variable  ( 42 )  is 
appropriate.  YHiile  nonlinear  scalings  have  been  encountered  before  in 
asymptotics,  even  in  combustion  (Sec.  II. 5),  an  exponential  one  appears  to 
be  novel.  It  is  fundamental  to  the  process  considered  here.  In  view  of  its 
role  in  thermal  runaway  it  is  not  surprising  that  p  is  the  other  appropria 
variable  for  the  explosion  phase.  In  terms  of  a,n  the  governing  equation 
becomes 

(63)  9T/3a  +  6[(n/2)3T/9n  -  32T/3p2]  =  (T2/Y0)(K0-T)exp(G/T0  -e/T  -0a). 


the  0(l)  values  of  0  and  n  for  which  the  equation  is  valid  correspond 
to  an  exponentially  thin  region,  the  hot  spot,  that  rapidly  gets  thinner 
as  time  increases.  The  focussing  as  thermal  runaway  is  approached  therefore 
continues  during  the  explosion  stage. 

If  the  reaction  term  in  equation  (63)  is  to  play  a  role,  we  must  write 

(65)  T  =  TQ/(1  -  cTq)  +  G"1T1(n,a)  +  o(e~'L) 

[cf.  the  result  (43)].  The  variable  a  only  appears  as  a  parameter  in 
the  equation  for  so  that  its  solution  may  be  written 
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(66)  T1  =  Tq{^(ti)  -  *n[(l  -  T0a)(Y0  -  H0T0a)/Y0] >/(l  -  TQa)2 
where 

(67)  -  ( ri/2)^/*  +  e^  =  1  with  i|)'(0)  =  0. 

This  is  identical  to  the  problem  (6b)  for  <t>Q>  indeed,  matching  the 
expansions  (36)  and  (65)  under  the  transformation  (42)  demands 

(68)  tj)  5  <j>Q. 

It  follows  that  the  boundary  layer  at  x  =  0  in  the  explosion  phase  is 
completely  determined  by  the  focussing  effect  in  the  induction  phase  and 
indeed  inherits  the  spatial  structure  of  that  effect. 

The  explosion  is  so  rapid  that  the  temperature  in  the  slab  away  from 
x  =  0  stays  sensibly  frozen  at  the  value 

(69)  T  =  TQ  +  0_1TQ(-2Xnx  -  in$  +  A)  +  o(e_1) 

attained  during  the  induction  phase.  Clearly  this  does  not  match  the 
expansion  (65),  even  to  leading  order.  An  intermediate  expansion  is  needed 
to  describe  the  structure  left  behind  by  the  rapidly  shrinking  hot  spot. 

The  breakdown  of  the  approximation  (69)  when  in  x  =  0(0)  suggests  that 

(70)  X  =  -204  in  x  =  0-  0"1  in(n2/S) 

is  an  appropriate  variable  to  use  with  a.  In  fact,  the  intermediate 
expansion  is  independent  of  0  ,  which  means  that  the  residual  structure 
is  stationary  too.  Details  have  been  given  by  Kapila  (1980);  we  merely 
note  here  that  the  leading  term  TQ/(l  -  T^X)  ranges  between  those  in 
the  expansions  (65)  and  (69)  as  X  decreases  from  0  to  0. 


Within  the  hot  spot  the  expansion  (65)  is  valid  until  0  approaches 
Yq/KqTq,  i.e.  T  approaches  .  If  a  gaseous  reactant  is  being  ecnsidere 
and  II  differs  appreciably  from  T^,  a$  it  often  does  in  practice,  the 
fluid  mechanics  can  no  longer  be  ignored  through  adoption  of  the  constant- 
density  approximation  (Kassoy  &  Poland  19Sob).  Otherwise  one  continues  with 
Kapila  (1980)  to  subsequent  phase  where  the  time  variable  (1»5)  is  appropriat 
that  is  not  suggested  the  expansion  (65)  as  far  as  we  have'  taken  it 
but  by  the  next  term.  The  corresponding  space  variable  is 


(71) 


t  -  S1/2  */e1/2 


since  that  ensures  a  diffusion  term  in  the  equation 
(72)  S(j>/3fr  -  3<j«/35  =  -$e<'1 


governing  the  coefficient  T.,(5,tj)  in  the  generalization  of  the  expansion 
(Ml).  The  initial  condition 


(73)  4>  =  -  in(-w)  -  Anin(-w)  +  ^(p)  +  o(l)  as  -»  with 

1  /2 

n  =  £/(-w )  fixed 

comes  from  matching  with  the  explosion  phase  while  the  boundary  condition 

2  2 

(7*0  <J>  -  -  in  5  -  inin  5  +  A  +  o(l)  as  g  ■>  «>  with  v  fixed 


comes  from  matching  with  the  intermediate  expansion,  neither  of  which 
will  be  derived  here;  in  addition 


(75)  =  0  for  5=0. 


These  conditions  fee***  the  origin  of  w 


undetermined,  as  in  the  homogeneous 


problem,  but  otherwise  ensure  a  unique  solution. 


j 

J&3 


■— »p -  f-;> -*  *■ r7r^r^^4*?i'‘w^  ”**  "i 
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Fig.  3  shows  the  picture  which  emerges  from  the  numerical  solution 
(Kapila  1980).  Since  the  initial  condition  (73)  is  effectively  independent 
of  5  when  that  variable  is  0(l),  the  curves  at  first  have  a  plateau  on 
the  left  which  is  continually  eroded  by  the  boundary  condition  (7^).  There 
comes  a  time,  however,  when  the  erosion  is  stopped  and  the  plateau  begins 
to  reestablish  itself  over  ever-increasing  distances,  i.e.  an  incipient 
deflagration  wave  is  formed.  The  focussing  effect  is  then  being  limited 
by  reactant  depletion  and  the  hot  spot,  now  fast  approaching  the  temperature 
Hq,  begins  to  spread  as  the  reaction  zone  detaches  itself  from  the  center 
and  starts  to  propagate  across  the  cold  unburnt  slab,  leaving  behind 
almost  depleted  reactant.  Most  of  what  was  the  final  relaxation  phase  in 
the  homogeneous  problem  is  now  taken  up  by  this  transition  from  hot  spot 
to  flash-through. 

The  remainder  and  more  is  swallowed  up  by  the  propagation  phase,  which 

takes  much  longer  than  the  transition  phase  but  is  still  very  brief.  Indeed 

-1  0/TO 

the  result  (XI.22)  shows  that,  for  an  0(0  e  )  value  of  D,  the 
deflagration  wave  takes  a  time  0(e€^^)  to  reach  the  surface  of  the  slab. 
This  suggests  using  the  variables 


(76) 


s  «  (x#(s)  -  xl/ee^2,  s  =  (T  -  x0)/e€‘ 


,1/2 


to  investigate  the  flash-through,  where  x,  is  measured  from  the  location 
x#  of  the  flame  front.  The  fact  that  0(9  changes  in  5  measure  the 
same  distances  as  0(l)  changes  in  5,  i.e.  that  the  flame  structure  will 
be  described  on  the  scale  of  5,  supports  this  suggestion. 

To  trace  the  progress  of  the  deflagration  wave  through  the  intermediate 
region  into  the  main  part  of  the  slab,  note  that  the  region  is  both 


,  ,  ,  r,T  -mv-iii 


-2k- 


stationary  and.  much  thicker  than  the  preheat  zone  of  the  wave.  Far  ahead 
of  the  flame  sheet  (on  the  scale  of  the  preheat  zone)  at  the  instant  it 
reaches  the  position  X  the  temperature  is  TQ/ ( X  -  TQX)  while  everywhere 
"behind  the  temperature  is  According  to  the  result  (il.22)  the  flame 

velociy  is  therefore  given  by 

(77)  x*(s)  =  H*(l  -  TQX)/C tQ  -  HT0X). 

As  X  decreases  from  Y^/K^T^  to  zero  through  the  intermediate  region, 

the  wave  decelerates  rapidly  from  infinity  to  v26  H^/Y '  ,  the  velocity 

at  which  it  flashes  through  the  main  part  of  the  slab.  (The  flame  moves 

more  rapidly  the  closer  the  unburnt  mixture  is  to  the  flame  temperature 

Hq  because  the  heat  of  combustion  is  able  to  preheat  a  larger  amount.) 

The  propagation  phase  ends  when  the  preheat  zone  reaches  the  surface. 

In  the  ensuring  relaxation  phase  the  wave  rapidly  decelerates  to  zero, 

thereby  establishing  a  steady  state  with  the  reaction  zone  just  inside  the 

slab.  Details  are  given  by  Kapila  (l98o);  we  merely  note  that  the  relaxation 

3/2 

take  a  very  short  time  0(06'  )  to  be  completed  [of.  the  0(£)  relaxation 
time  for  the  homogeneous  explosion]. 

We  do  not  go  into  details  of  the  relaxation  phase  because  it  is  a 
creature  of  the  unphysical  boundary  condition  on  Y.  If  instead  of  maintaining 
the  reactant  at  the  surface  iff  is  prevented  from  diffusing  across  the  surface, 
then  nothing  changes  until  the  relaxation  phase  since  Y  +  T  is  constant  to 
leading  order  everywhere.  At  that  time  the  wave  does  not  come  to  rest, 
but  is  quenched  by  the  cold  surface  (Buckraaster  1979) >  leaving  behind 
reactant-depleted  material  at  temperature  Hq  except  for  a  small  amount  of 
unburnt  reactant  near  the  surface.  Since  the  boundary  continues  to  be 
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maintained  at  Tq,  there  follows  a  relaxation  phase  of  duration  0(l) 
during  which  essentially  pure  heat  conduction  loweis  the  temperature 
everywhere  to  T^  again  and  the  residual  reactant  spreads  across  the 
slab,  being  slowly  burnt  in  the  process. 

Our  analysis  describes  the  formation  of  a  hot  spot  from  which  emerges 
a  deflagration  wave  that  rapidly  consumes  the  combustible  material. 

These  are  chai-acteristics  of  actual  explosions,  so  that  the  model  has 
clearly  retained  the  essential  physics.  For  a  solid,  modifications  are 
necessary  to  account  for  infinite  Lewis  number  but  these  only  change  details. 
For  a  gas,  the  neglect  of  the  fluid  mechanics  has  undoubtedly  excluded 
significant  effect  (Kassoy  &  Poland  1980b),  in  particular  the  evolution  of 
the  deflagration  into  a  detonation.  The  real  value  of  the  discussion, 

‘however,  lies  in  its  elucidation  of  the  extremely  rapid  transients  typical 
of  the  explosion  process,  one  of  the  most  challenging  aspects  of  combustion 
theory. 

5. '  Ignition. 

Thermal  explosion  is  triggered  in  a  combustible  material  by  the  build-up 
of  heat  which  cannot  escape  in  sufficient  amounts  through  the  boundary. 

By  contrast,  ignition  is  caused  by  some  external  source  of  energy.  For 
example,  part  or  all  of  the  boundary  may  slowiy  be  raised  to  a  high 
temperature  so  that  reaction  starts  nearby.  Such  a  problem  was  considered 
in  Sec.  2  for  a  slab.  (Sudden  elevation  of  temperature  is  of  more  interest.) 
Alternatively,  a  radiative  heat  flux  may  be  applied,  a  mechanism  that  is  some¬ 
times  a  factor  in  the  spread  of  fires. 


•x.'sSii  i  v. 
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Thermal  runaway  is  not  only  a  characteristic  of  explosion  but  also 
of  ignition.  As  is  clear  from  Sec,  4,  if  the  goal  is  limited  to  calculating 
the  time  to  runaway,  then  only  the  induction  stage  need  be  considered. 
Indeed,  that  is  as  far  as  the  analysis  has  been  taken  by  anyone  to  date. 

Our  intention  is  not  only  to  draw  attention  to  ignition  problems  and 

what  needs  to  be  done  (problems  which  only  make  sense  in  the  context  of 

« 

large  activation  energy)  but  also  to  disclose  another  way  in  which  runaway 
can  occur. 

The  specific  problem  we  shall  examine  is  that  of  ignition  of  a  half 
space  by  a  constant  energy  flux  (Linan  &  Williams  1971).  Modifications  for 
a  slab  with  its  other  surface  insulated  (say)  are  minor  and  will  be  obvious. 
Since  the  discussion  does  not  go  beyond  thermal  runaway,  reactant  depletion 
may  be  neglected;  so  that  only  the  temperature  equation  (17)  need  be 
considered,  i.e. 

(78)  3T/3t  -  32T/3x2  =  De-0^1’ 

in  the  present  geometry.  The  boundary  and  initial  conditions  are 

(79)  3T/3x  =  Tq  >  0  for  x  =  0,  T  =  Tq  for  t  =  0  and  as  x  -» 

if  the  half-space  lies  in  -«  <  x  0.  Initially  the  temperature  is  uniform 
so  that,  in  the  absence  of  external  stimulus  (e.g.  when  the  surface  is 


insulated) ,  there  will  eventually  be  an  explosion  if  D  is  large  enough. 
However,  we  are  concerned  with  the  (earlier)  occurrence  of  runaway  due  to 
the  constant  heat  flux  applied  at  the  boundary. 
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To  establish  some  preliminary  ideas  define  T  ,  a  temperature 
characterizing  the  reactivity  of  the  material,  by 

(80)  D  =  A(e)ee^Tr 

where  A(0)  is  algebraic  in  0.  Clearly  Ty  should  be  greater  than  Tq» , 

otherwise  there  will  be  reaction  from  the  start  (0  <  T  <  T.)  or  no  reaction 

r—  0 

at  all  (T  <  0).  The  equation  governing  T  is  then 

(81)  3T/8t  -  32T/8x2  =  A(@)exp(0/Tr  -  e/T) 

where,  at  times  sufficiently  small  for  T  to  be  less  than  T  .  everywhere, 
the  right-hand  side  is  negligible.  The  solution 

(82)  T  =  Tq  +  T^[2/t/ir  exp(-x2At)  +  x  erfc  (|x[/2/t)] 

satisfying  the  conditions  (79)  increases  monotonically  in  t  but  decreases 
monotonically  in  |x|,  so  that  the  value  Tf  is  first  attained  at  x  =  0, 
the  time  being 

(83)  tr  =  7f(Tr  -  Tq)2At^2. 

At  that  instant  the  reaction  term  must  be  reinstated.  [To  obtain  t  in 

r 

terms  of  D,  correct  to  leading  order,  insert  Tr  =  Q/i n  D. ] 

An  accurate  estimate  of  the  runaway  time  must  involve  A(@),  i.e.  the 
reaction  (as  in  Sec.  *0.  Nevertheless,  the  notion  that  the  temperature 
rises  purely  by  heat  conduction  until  it  reaches  some  specific  value 
determined  by  the  magnitude  of  D,  whereupon  reaction  becomes  important 
enough  to  produce  runaway  within  a  time  0(0  ^),  is  central  to  the 


development. 


(90) 
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the  temperature  being  too  far  from  T for  there  to  be  any  reaction  (to 
any  order),  Matchings with  the  inner  and  outer  expansions  and  with  the 
solution  (82)  require 


(97) 

Tx  =  0,  3TX/3X  =  Tj 

for 

0 

11 

X 

(98) 

T  =  TJX  +  o(l)  as 

X  -v 

-00  and  as  a  ->-  -». 

There  is  one  boundary  condition  too  many  but  that  does  not  prevent  there 
being  a  solution 


(99) 


*i  *  V- 


Continuing  to  T„  yields 

(100)  (d/ 3a  -  32/3X^)T2  =  0, 

(101)  T2  =  T1q,  3T2/3X  -  (CT^/T^)exp(T10/T^)  for  X  =  0, 

(102)  T2  =  T^X2  +  2a)/2^T  +  o(l)  as  X  -*■  -»  and  as  a 

unlike  that  for  T^,  this  problem  is  not  overdetermined  since  T^q  has  yet 
to  be  found.  Indeed  its  elimination  from  the  boundary  conditions  (101)  gives 
the  single  nonlinear  condition 


(103)  3T2/3X  =  (CT2/Tj)exp(T2/T2)  for  X  =  0. 


/ill  parameters  can  be  purged  by  writing 

r 


( 10U ) 


i  T2  =  TjfX2  +  2c)/2\'rit^'  +  T2  %,  X=  (w^/t'2)1^?, 
a  =  (irtrT^2)l/2[a  +  £■  ^(t^M^C4)] 


'}”.****■£  V 
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to  obtain 

(105) 

(106) 
(107) 


O /8a  -  92/s3?)t1  =  0, 


ST^/SX  =  exp(cr  +  T9)  for  X  =  0, 


Tr,  =  o(l)  as  X  ■*■  -»  and  es  5  +  -« 


Numerical  solution  of  this  problem  shows  that  Tg  becomes  unbounded  at  a 
definite  time 


(108) 


a  =  -.1*31 


(Lilian  &  Williams  1971).  In  terms  of  the  original  time  variable,  runaway 
therefore  occurs  at 

» 

(109)  tr  +  6“1(TitrT^/T'2)1/2[-.431  +  £  to  (T^/irt^C11)]. 


The  mathematical  problem  for  runaway  can  therefore  assume  different 

forms.  The  nonlinear  driving  term  is  in  the  governing  equation  (55)  and 

acts  on  0(l)  time  and  length  scales.  Here  it  is  in  the  boundary  condition 

(106)  and  acts  over  a  time  0(9  )  at  distances  0(6”  ).  Moreover  the 

occurrence  of  runaway  depends  on  the  reactivity  of  the  material  in  some 

cases  but  not  in  others.  The  explosion  of  Sec.  4  requires  Ty  =  Tq  and 

materials  of  lower  reactivity  (Tr  >  TQ)  will  never  explode;  on  the  other 

hand,  the  ignition  considered  here  will  occur  for  any  Tf  >  Tq,  as  will  the 

explosion  of  Sec.  3  (albeit  after  an  exponentially  long  warm-up  period)  if 

T  is  also  less  than  H-. 
r  0 


\ 
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6.  Explosion  of  Separated  Reactants. 

If  two  parallel  streams,  one  a  combustible  mixture  and  the  other 
a  hot  inert,  are  allowed  to'  come  into  contact  then  eventually  the  mixture  will 
be  ignited  by  the  inert.  Such  is  the  steady  problem  considered  by  Marble 
&  Adamson  (195*0*  Here  we  shall  treat  a  somewhat  different  problem  having 
some  of  the  same  features  but  closer  to  the  thermal  explosion  investigated 
•  in  Sec.  4.  Mathematically  it  is  identical  to  a  special  case  of  the  unsteady 
mixing  and  explosion  of  two  initially  separated  reactants  considered  by 
Linan  &  Crespo  (1972).  Indeed  we  shall  retain  the  notation  of  the  unsteady 
problem  so  as  to  use  results  from  Sec.  h  without  translation.  One  feature 
missing  from  the  Mar b 1 e -Adam s on  problem  is  the  diffusion  flame  which  manifests 
itself  as  a  Burke-Schuman  flame  sheet  sufficiently  far  downstream. 

Consider  a  semi-infinite  flat  plate  (Fig.  *0  separating  fast  parallel 
streams  containing  oxidant  X  and  fuel  Y  that  come  into  contact  at  the 
trailing  edge  and  then  interdiffo.se  to  form  a  combustible  mixture  down¬ 
stream.  For  simplicity  the  speeds^ reactant  concentrations  and  temperatures 
of  the  streams  will  be  taken  equal.  (Linan  &  Crespo's  problem  is  equivalent 
to  one  allowing  for  unequal  concentrations  and  temperatures.)  For  small 
values  of  t  no  significant  reaction  occurs;  for  large  values  there  is 
chemical  equilibrium,  with  a  Burke -Schumann  flame  sheet  separating  a 
region  in  which  X  is  zero  from  one  in  which  Y  is  zero.  How  the  com¬ 
bustion  field  evolves  from  the  initial  regime  to  the  final  one  as  t 
increases  is  the  main  question.  Such  a  transition  arises  in  many  combustion 
processes,  another  example  being  the  unsteady  ignition  of  a  fuel  drop,  so  that 
a  m&l  .ematical  description  is  of  considerable  interest.  Moreover, 
experiments  have  been  carried  out  under  conditions  similar  to  those  in 
Fig.  *4  [Liebman,  Corry  &  Perloe  (1970)). 


The  essential  details  of  the  combustion  field  can  be  uncovered  by 
activation-energy  asymptotics,  as  was  shown  by  Lilian  &  Crespo  (1972).  Their 
discussion  is  highly  intuitive  since  asymptotic  methods  were  not  well 
developed  at  that  date,  but  nevertheless  they  accurately  described  the 
main  features  of  their  more  general  problem.  The  present  discussion  is 
more  systematic  and  the  description  richer;  only  for  the  sake  of  simplicity 
is  the  problem  specialized. 

The  high  velocity  of  the  streams  enables  longitudinal  diffusion  to 
be  neglected  (cf.  Sec.  VIII. 6);  in  addition  we  adopt  the  constant-density 
approximation.  Taking  unit  Lewis  numbers  then  leads  to  the  governing 
equations 

(110)  (3/3t  -  32/3x2)T  =  -2(3/3t  -  32/3x2)X  =  -2(3/3t  -  32/3x2)Y  =  D(XY/Y2)e 

» 

if  the  reaction  is  first-order  in  both  fuel  and  oxidant.  Such  equations, 
being  parabolic,  require  only  initial  condition;  we  shall  take 

.rx  ,  ro  for  x  <  0 

(111)  T=T,X=SUY=S  at  t  =  0 

,  ^YQ  for  x  >  0 

with 


(112) 


xo  =  yo 


to  guarantee  symmetry  about  the  t-axis,  a  simplification  that  does  not 
exclude  any  essential  feature.  The  variable  t  is  time-like,  so  that 
corresponding  terminology  will  be  adopted;  when  t  is  interpreted  as  time, 
equations  (110)  govern  Lifi£n  &  Crespo's  one-dimensional  unsteady  problem 
of  two  half-spaces  of  fuel  and  oxidant  brought  suddenly  into  contact  at 
the  initial  instant.  (Their  generalization  consists  of  different 
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temperatures  for  x  ^  0  and  XQ  ^  YQ;  they  also  take  variable  density 
into  account  by  means  of  the  Howarth-Dorodnytsin  transformation ) 

The  linear  combinations 

(113)  G  =  T  +  2X,  H  =  T  +  2Y 

satisfy  pure  heat  equations  everywhere  with  step  functions  for  initial 
values;  whence  follow  the  Shvab-Zeldovich  relations 

(lilt)  G  «  Hq  -  Yq  erf  n  ,  H  =  HQ  +  YQ  erf  n  (Hq  =  Tq  +  YQ) , 

where  the  similarity  variable 
(115)  n  =  x/2/t 

arises  naturally.  We  now  need  only  consider  the  temperature  equation 
which,  in  terms  of  t  and  n  as  independent  variables,  becomes 

f  .V  3T  lp  3T  1  32T  _  +T.fYY/v2s  -6/T 

(H6)  t  v  -  tD(XY/Y0)e 

where  X  and  Y  have  just  been  determined  as  functions  of  T  and  n  . 

Behavior  for  small  and  large  values  of  t  can  be  inferred  from  this  form 

of  the  equation,  for  any  value  of  0  ,  by  noting  that  the  effective 
•  * 

Damkohler  number  is  tD:  for  t  small  the  combustion  is  nearly  frozen 
whereas  for  t  large  it  is  near  equilibrium. 

In  the  limit  t  ->  0  we  find 

(117)  T  =  T0,  X  =  YQ(1  -  erf  n)/2,  Y  =  YQ(l  +  erf  n)/2. 

At  any  finite  value  of  x  when  b  is  small  there  is  a  combustible 
mixture,  formed  by  interdiffusion  of  the  reactants,  in  which  no 
significant  reaction  has  taken  place.  We  may  therefore  expect  a  thermal 


explosion  similar  in  many  respects  to  the  spatially  inhomogeneous  one 
discussed  in  Sec.  4. 

In  the  limit  t  •*  «  we  must  have 


(118) 


XY  =  0 


and  then  the  relations  (ll4)  show  .that 


(119)  T  = 


Hq  +  Y0  erf  n,  j  ~Vrf  * 


*0  *  Yo  erf 


0 

.  x=l 

n,  ^  o 


,  Y  =i 


0  for  x  <  0, 


Y0  erf  n  for  x  >  0. 


There  is  a  Burke-Schumann  flame  sheet,  corresponding  to  infinite  Pamkohler 
number  (cf.  Sec.  VI. 2),  located  on  the  t-axis  at  a  temperature 


(120) 


T*  =  no¬ 


on  the  scale  cf  n  the  flame  sheet  is  a  discontinuity  so  that,  to 
investigate  its  structure,  we  introduce  the  variables 


(121)  tl/3n  =  x/2tl/6,  ©=  /TTtl/3(T«  -  T)/2Yq 


and  seek  a  solution  for  which  is  a  function  of  S  alone.  This 
choice  of  variables  ensures  a  meaningful  balance  of  diffusion  and  reaction, 
i.e.  a  distinguished  limit  of  equation  (ll6)  in  asymptotic  parlance.  The 


limit  is 


(122) 


d2<57dc2  =  C(G2  -  ?2)  with  C  =  2De  ^  Ql(^\ 


and  matching  with  the  outer  solutions '( 119)  requires  the  boundary  conditions 


(123) 


d0/d?  =  ±1  +  oil)  as’  ±oo  . 
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The  problem  (122,123)  is  that  for  the  classical  Burke-Schumann  structure 
(cf.  Sec.  VI. 2). 

The  evolution  from  the  frozen  regime  to  the  equilibrium^  regime 
turns  out  to  be  remarkably  simple  in  the  limit  8  -*■  ».  Thermal  runaway 
occurs  in  the  diffusion-generated  mixture  at  a  point  on  the  line  of 
symmetry  and  a  hot  spot  develops  just  as  for  the  spatially  inhomogeneous 
explosion  of  Sec.  4.  As  the  temperature  in  the  hot  spot  approaches  its 
ultimate  value  (corresponding  to  complete  consumption  of  the  reactants) 

a  symmetrical  pair  of  deflagration  waves  is  formed,  which  penetrate  the 
large  temperature  gradients  left  behind  by  the  focussing  hot  spot  and 
flash  through  the  layers  of  frozen  mixture  consuming  all  of  the  deficient 
reactant.  All  this  happens  within  a  range  of  t  (including  the  runaway 
point)  that  is  exponentially  smaller  than  the  scale  of  mixing. 

A  detailed  description  starts  with  the  induction  phase,  during 
which  the  temperature  rises  from  the  initial  value  TQ  by  only  an 
0(0  ■*■)  amount.  The  corresponding  small-disturbance  equation 

(124)  94-/ 3t  -  3 24>/3x2  =  6(1  -  er^rDe* 

differs  from  that  (55)  for  an  initially  uniform  slab  only  in  having  an 

extra  factor  on  the  righthand  side;  if  we  now  restrict  the  discussion 

A 

to  the  region  x  >  0,  the  boundary  conditions 

(125)  3$/3x  =  0  for  x  =  0,  4=  0  for  t  -  0 

are  the  same.  As  there  the  solution  becomes  unbounded  at  a  point 

(126)  t  =  t  on  x  =  0 

e 

which  can  be  determined  numerically  (Linan  ?<  Crespo,  1972). 
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The  behavior  of  $  as  t  +  t  e  is  given  by  equation  (59)  since  n 

is  0(t  -  t)1^2  in  the  variables  t  and  n  •  The  function  <j>A  is 
c  u 

changed  because  the  constant  A  in  its  asymptotic  expansion  (6l)  is  different, 
being  determined  by  the  asymptotic  form  (62)  of  the  solution  to  a  different 
problem.  (The  value  of  A  has  never  been  computed.) 

The  explosion  phase  is  governed  by  the  equation  (63)  with  right-hand 

* 

side  replaced  by 

(127)  (T2AY2)(Ho-T)2exp(0/To  -8/T  -6a) 

<• 

since  n  is  effectively  zero  there.  The  reactant  concentrations  are 
changed  from  Y^/ 2  to  zero  as  T  increases  from  T^  to  H^.  The 
expansion  (65)  is  still  valid,  but  now 

* 

(128)  Ta  =  T24(t,)  -  2*n[(Y0  -  H0T0a)/2Y0l}/(l  -  TQo)2 

where  again  \Ji  satisfies  the  equation  (67)  and  hence  the  identity  (68) 
with  the  lew  <}>q.  Between  the  hot  spot  described  by  this  expansion  and 
the  essentially  frozen  combustion  field  on  either  side  there  is  a 
stationary  intermediate  structure  on  the  scale  (70),  the  leading  term 
of  which  is  again  TQ/(l  -  TQX). 

The  transition  phase  which  follows  is  governed  by  the  modification 

(129)  8<*>/8w~  32<j>/9§2  =  -<j>2e^ 

of  the  equation  (72),  provided  the  definition  {b6)  is  replaced  by 

(130)  f  =  (4Y202/H2T2)exp(-9Yo/KoTo). 

The  modification  arises  because  the  hot  spot  contains  the  reactants  in 
stoichiometric  proportion,  which  results  in  a  truly  second-order  reaction 


rather  than  an  essentially  first-order  one.  The  initial  condition  (73)  is 
replaced  hy 

(131)  <j>  =  -£n(~to)-2$,n£n(-ii))  +  i|i(n)  +  oil)  as  u  +  -« 

with  n  =  C/(“u>)1//2  fixed 

.  while  the  boundary  condition  (7*0  is  replaced  by 

p  2 

(132)  <}>  =  -£n£  -  2£ntn  £  +  A  +  0(l)  as  5  +  »  with  u  fixed. 

It  appears  that  no  numerical  investigation  of  the  problem  (129,  131,  132) 
has  yet  been  undertaken,  but  similarity  to  the  problem  (72,73,7*0  leaves 
little  doubt  that  it  would  yield  a  transition  from  a  focussing  hot  spot 
to  an  incipient  deflagration  wave. 

The  wave  must  first  traverse  the  intermediate  structure  left  by  the 
shrinking  hot  spot.  The  situation  is  the  same  as  that  in  Sec.  *1  except 
that  ahead  there  is  a -mixture  of  fuel  and  oxidant  in  stoichiometric 
proportion,  i.e.  X  =  Y  [cfi  equations  ( llU )  with  n  =  0] ,-  instead  of 
a  single  reactant.  It  follows  that  the  variables  (76)  are  again 
appropriate,  if  the  new  definition  (130)  of  6  is  used,  but  that  the 
velocity  (77  bis)  must  be  replaced  by 

H^t  1  —  T  p( 

(133)  W  «/fi— exp(Y  0/2II  T  )(— ~-^) 

Y0e2  °  °  °  Wo* 

in  accordance  with  the  result  (II.52).  As  X  decreases  from  Yq/HqTq 
to  zero  through  the  intermediate  region  the  wave  decelerates  rapidly  from 


infinity  to 


the  velocity  with  which  it  enters  the  near-frozen  external  field. 


For  the  slab  considered  in  Sec.  k,  the  wave  then  continues  at  the 
same  speed  to  the  surface,  but  here  it  encounters  non-uniformities  in 
available  enthalpy  which  produces  0(l)  changes  in  its  temperature  and  hence 
exponentially  large  changes  in  its  velocity.  These  non-uniformities  have 
been  formed  by  interdiffusion  of  the  reactants  up  to  t  =  Te/^>  i*e. 

(135)  T  s  Tq,  X  =  YQ(1  -  erf  ne)/2,  Y  =  YQ(l  +  erf  n^/2 
where 

(136)  ne  =  x/2/te, 

and  have  no  time  to  change  during  transit  of  the  wave. 

The  situation  is  similar  to  that  in  the  intermediate  region:  the 
combustion  field  is  stationary  and  has  a  much  larger  scale  than  that  of 
the  preheat  zone  of  the  flame.  Far  ahead  of  the  wave  on  the  latter  scale 
the  temperature  and  mass  fractions  have  the  distributions  ( 135 )  at  the 
instant  the  wave  reaches  the  position  x;  while  behind 

(137)  T  =  HQ  -  YQ  erf  ng,  X  =  0,  Y  =  YQ  erf  ne  for  x  >  0. 

To  determine  the  velocity  of  the  wave  we  use  the  result  (ll.5l)  which,  as 
was  pointed  out  there,  holds  both  near  and  far  from  stoichiometry.  Thus 

(138)  W  *  /6(To/YQ02)(l-erffie)";L(Ko-Yoerfne)[(Ho-Yoerfne)2+0Yoerfne]l/2eaC 
where 

(139)  a  =  Yq(1  -  erf  ^)/2T0(H0-Y0  erf  ng) 
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is  the  wave  speed,  which  clearly  reduces  to  the -value  (134)  for  ng  =  0. 
The  result  cannot  be  valid  as  n  -*■  «  since  the  propagation  is  then 
associated  with  significant  changes  in  t;  well  before  then,  however,  the 
formula  (ll.5l)  fails  because  Jg  (i.e.  the  oxidant  flux)  becomes  small. 
The  correct  result  has  not  yet  been  obtained,  perhaps  because  it  would 
add  very  little  to  the  picture.  Note  that  equilibrium  with  a  Burke- 
Schumann  flame  sheet  holds  at  any  finite  x  immediately  after  the  explo¬ 
sion:  the  combustion  field  (137)  is  precisely  of  the  ultimate  form  (119b). 

The  picutre  which  has  emerged  is  remarkably  simple,  consisting  of 
regions  of  frozen  chemistry  and  equilibrium  with  reaction  zones  between 
them,  the  hallmark  of  activation-energy  asymptotics.  To  secure  the 
essential  features  very  little  analysis  was  required  beyond  that  for 
the  thermal  explosion  of  Sec.  4. 
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